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Abstract 

We consider the possibility that the dark matter is coupled through its mass to a scalar 
field associated with the dark energy of the Universe. In order for such a field to play a role 
at the present cosmological distances, it must be effectively massless at galactic length scales. 
We discuss the effect of the field on the distribution of dark matter in galaxy halos. We show 
that the profile of the distribution outside the galaxy core remains largely unaffected and 
the approximately flat rotation curves persist. The dispersion of the dark matter velocity 
is enhanced by a potentially large factor relative to the case of zero coupling between dark 
energy and dark matter. The counting rates in terrestrial dark matter detectors are similarly 
enhanced. Existing bounds on the properties of dark matter candidates can be extended to 
the coupled case, by taking into account the enhancement factor. 

PACS numbers: 95.35.+d, 98.35. Gi, 98.80.Cq 



1 



According to our present understanding of the matter content of the Universe, its energy 
density is dominated by dark energy and dark matter. The usual assumption is that the dark 
energy is associated with the potential of an evolving scalar field [1, 2], while the dark matter 
is composed of weakly interacting massive particles. Usually these two sectors arc assumed 
to interact only gravitationally. Attempts to resolve the coincidence problem (the comparable 
present contributions to the total energy density from the two components) are based on the 
presence of an additional interaction between the dark matter and the scalar field [3]. The 
cosmological evolution now depends on the potential of the field, as well as on the type and 
strength of the interaction [4]. For a scalar field to be relevant for the cosmological evolution 
today, so as to resolve the coincidence problem, the scale for its effective mass must be set by 
the present value of the Hubble parameter [3] . This means that the field is effectively massless 
at galactic length scales. Its coupling to the dark matter particles results in a long-range force 
that can affect the details of structure formation [5, 6]. 

The presence of a long-range scalar interaction in addition to gravity can modify the dis- 
tribution of dark matter in galaxy halos. In particular, it is not certain that the dark matter 
distribution can maintain a form that results in approximately flat rotation curves for objects 
orbiting galaxies. Recent simulations of structure formation [6] are not detailed enough at the 
scales of interest in order to resolve this question. 

We consider an interaction between the scalar field and the dark matter particles that can 
be modelled through a field-dependent particle mass. The action takes the form 

S = J d^xV^ [m^R - ^g^^'d^cj) d,(P - Ui(P)^ m{(t>{xi))dTi + S^, (1) 



where drj = J —g^i,{xi)dx^dx\ and the second integral is taken over particle trajectories. The 
contribution describes the standard baryonic component. Variation of the action with respect 
to (j) results in the equation of motion 



\/—g d^ d(f) 

where the energy-momentum tensor associated with the gas of particles is 

/ dn mi(j)(xi)) -4 

-9 



T^^ = ^ / dn m{cl>{xi)) ^^5('\x - Xi). (3) 



We are interested in static spherically symmetric configurations, with the scalar field varying 
slowly with the radial distance r. Our treatment is relevant up to a distance ri ~ 100 kpc beyond 
which the dark matter becomes very dilute. For r ^ ri we expect that (p quickly becomes 
constant with a value close to <i){ri) = (f)\. This is the value that drives the present cosmological 
expansion. Here we assume that the cosmological evolution of is negligible for the time scales 
of interest, so that the asymptotic configuration is static to a good approximation. 

We approximate: m{4>) ~ m{(po) + [dm{4>o) / d(j)] 6<p = rno + JTIqJ^, with cpQ the value of the 
field at the center of the galaxy (r = 0). We work within the leading order in Scf) and assume 
that m'/m ~ mg/mo for all r. Also dU/d(j) can be approximated by a constant between r = 
and r = GO. For the scalar field to provide a resolution of the coincidence problem, the two 
terms in the r.h.s. of eq. (2) must be of similar magnitude in the cosmological solution. This 
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means that dU/dcp must be comparable to (mo/mo)/3oo- We expect poo to be a fraction of the 
critical density, i.e. ^ 3 keV/cm^. On the other hand, the energy density in the central 
region of the static solution (r ^ 100 kpc) is that of the galaxy halo (~ 0.4 GeV/cm^ for our 
neighborhood of the Milky Way). This makes dU/d(f) negligible in the r.h.s. of eq. (2) for a 
static configuration. The potential is expected to become important only for r oo, where 
the static solution must be replaced by the cosmological one. Similar arguments indicate that 
we can neglect U relative to p. Also the scalar field must be effectively massless at the galactic 
scale. For these reasons we expect that the form of the potential plays a negligible role at the 
galactic level. Our analysis can be carried out with {7 = and is model independent. 

We treat the dark matter as a weakly interacting, dilute gas. We are motivated by the 
phenomenological success of the isothermal sphere [7] in describing the flat part of the rotation 
curves. We do not address the question of the density profile in the inner part of the galaxies 
(r ^ 5 kpc). We approximate the energy-momentum tensor of the dark matter as T'^ = 
diag{—p, p,p,p) with p(r) = p{r) (v^) = m(4>{r)) n{r) (v^). The dispersion of the dark matter 
velocity is assumed to be constant and small: (v^) ^ 1. The gravitational field is considered in 
the Newtonian approximation: qqq ~ 1 + 2$, with $ = O (mQS(p/mo). In the weak field limit 
and for p <^ p, the conservation of the energy-momentum tensor gives 



m'n 



p' = -p^'-p^{5<py, (4) 

mo 

with the prime on p, 5(p denoting a derivative with respect to r. Integration of this equation 
gives ~ no exp (— ^*/(v^) — {mQ/mo)S(j)/ (v^)) . In order to make the picture more complete we 
also allow for a pressureless baryonic component in the core. We assume that its energy density 
has the phenomenological profile pb{r) = ps f (f/fc), with f(x) a decreasing function of x. We 
do not discuss the physics that leads to such a profile. We include the baryonic contribution 
only in order to estimate its effect on the dark matter distribution in the outer regions of the 
galaxy. 

With the above assumptions we obtain the equations of motion 

^" + = 4i^ ^""^ " "'^'^^ + 4i^ ^ (^:) ' 

and 

i5<j>y' + - iScj)y = !^ po exp - aScP) , (6) 
r mo 

where M = {16irGi^)~^/^ is the reduced Planck mass, po = '"^o^o the energy density of dark 
matter at r = 0, a = l/(v^), and a = mo/(?n.o(v^)). We emphasize that, even though |$| <^ 1, 
the combination $/(v^), that appears in the exponent in the expression for the number density 
n, can be large. Similarly, the expansion of the mass around the value mo = m((^o) assumes 
the smallness of the dimensionless parameter |mo(5<^/mo|. However, the combination aScf) = 
(m,o5(/)/mo)/(v^), that appears in the exponent, can be large. 
A linear combination of eqs. (5), (6) gives 



d'^u 2du R r f z \ ^ 

where u = — a$ — a5(/>, z = l3r, Zc = fir^,, R = Pb/ Pq-, and = (1+k^)q;po/4M^. The parameter 
= 4M^ (mo/mo)^ determines the strength of the new interaction relative to gravity. For large 
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z all regular solutions approach the form expu = Ij . Another linear combination of eqs. (5), 
(6) gives 

d^v 2 dv K^R , f z\ 



dz"^ z dz 1 + \Zc/ 

with V = —K?a^ + a6(f>. For large z the solution of this equation is v = cq + ci/ z. Combining 
the two solutions gives the leading result for large r 



If the new interaction is universal for ordinary and dark matter, the experimental constraints 
impose <C 1. In this case, it is reasonable to expect a negligible effect in the distribution of 
matter in galaxy halos. However, if interacts only with dark matter, as we assume here, this 
bound can be relaxed significantly. 

According to cq. (9), a massive particle in orbit around the galaxy, at a large distance r from 
its center, has a velocity v^ = 2(v^) /(I + k^). We can use this expression in order to fix (v^) for a 
given value of k. The approximately flat rotation curves are a persistent feature even in the case 
that the dark matter is coTiplcd to a scalar field through its mass. The effect of the new scalar 
interaction is encoded in the factor k? . When this is small, the velocity of an object orbiting 
the galaxy is of the order of the square root of the dispersion of the dark matter velocity. If 
is large, the rotation velocity can become much smaller than the typical dark matter velocity. 

The Newtonian potential can be expressed as = —2 {u' + v') /v^. This allows us to write 
the rotation velocity as 

V \ 2 z f du dv\ , ^ 

vj =T = -2U + d^)- ^''^ 

We can also relate the parameters z and r through the expression 

^'^-\0.2MQ/pcy Ve kpc' ^ ^ 

In order to study the rotation velocity profile more quantitatively we integrate eqs. (7), (8) 
in the interval < z < 30. We normalize the Newtonian potential $ so that $(0) = 0. Also, 
according to our definition, 5(/>(ro) = 0. These relations imply u(0) = v{0) = 0. The initial 
conditions du(0)/dz = dv{0)/dz = guarantee that the solution is regular at r = 0. In fig. 1 
we present the solution of eqs. (7), (8) for various values of k and R. In all cases we use f{x) = 
exp(— x) and Zc = 0.5. We observe that the solutions do not differ substantially. Approximately 
flat rotation curves are expected for z ^ 5, consistently with observations [7]. This behaviour 
persists for other forms of the profile of baryonic matter, as long as dark matter dominates over 
baryonic matter for z ^ 3. The flat part of the rotation curves is largely insensitive to the 
structure of the core. The departure from spherical symmetry for the distribution of matter in 
the core of a spiral galaxy is also expected to be a minor effect. 

Independently of the value of k^, the interaction of dark matter with the scalar field asso- 
ciated with dark energy does not destroy the approximately flat profile of the rotation curves. 
Other considerations, however, could constrain the coupling = 4M^ (mo/mo)^. The dis- 
persion of the dark matter velocity is (v^) = (1 + k^)v^/2. For a value of Vc deduced from 
observations, (v^) increases with k. For sufficiently large k, it seems possible that Vc may exceed 
the escape velocity from the galaxy. It turns out, however, that this is not the case. Outside the 
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Figure 1: The shape of the rotation curves for various values of k and R. 



core of the galaxy and for r ^ ri, the binding potential for a dark matter particle is ^ + {di/a)5(j). 
For large r, eq. (8) implies that v = —K^a^ + a5<j) = constant. The binding potential becomes 
= (1+K^)v^ ln(r/ri), where we have omitted an overall constant. For a particle at a dis- 
tance from the center of the galaxy, the escape velocity becomes v^^^ = 2(l-|-K^)[ln(ri/r^,)-|-l]. 
The value of \esc is larger than the standard one [8] by a factor (1 + k?), so that (v^) remains 
substantially smaller than v^^^ for ^ ri. A particle that does not interact with the scalar 
field is bound only by the the potential ^. However, the scale of its velocity is set by Vc, so that 
again it cannot escape. 

Our discussion supports the conclusion that the presence in galaxies of a scalar field related 
to dark energy does not have any disastrous phenomenological consequences. On the contrary, 
this scenario may be much more interesting for dark matter searches than the conventional one. 
Depending on the value of k, the typical velocity of dark matter particles can exceed significantly 
the observed rotation velocity (~ 220 km/s for the Milky Way). As the density profile is not 
significantly modified, the estimated local energy density of dark matter remains the same as 
in the case with = 0. It is ~ 0.4 GeV/cm^ for our neighborhood of the Milky Way. As a 
result the fiux of dark matter particles towards a terrestrial detector is larger roughly by a factor 
(1 + K^)^/^ relative to the k? = Q case. A detailed study must take into account the motion 
of the Earth through our galaxy. However, the Earth velocity is of the order of Vc, and gives 
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only a modest correction. The cross section for the elastic scattering of halo particles by target 
nuclei through weak or strong interactions is largely independent of the particle velocity [9] . The 
leading effect of a non-zero value of k is that the counting rates, that are proportional to the 
velocity, arc increased by the factor (1 + k^Y^"^ . This makes the dark matter easier to detect. 
Existing bounds on dark matter properties from direct searches can be extended to include the 
case of non-zero k. The bound on the cross section for the interaction of dark matter with the 
material of the detector must be strengthened by (1 -|- k^)^/^. 

The allowed range of k is limited by the observable implications of the model that describes 
the dark sector. Strong constraints arise from the measured fluctuations of the cosmic microwave 
background. The dependence of the mass of dark matter particles on an evolving scalar field 
during the cosmological evolution since the decoupling is reflected in the microwave background. 
The magnitude of the effect is strongly model dependent. In the models of rcf. [4, 10], which 
assume an exponential dependence of the dark matter mass on the field, the observations result in 
the constraint ^ 0.01. In the model of refs. [5, 6] the dark matter mass is a linear function of 
the field. The induced scalar interaction among dark matter particles is screened by an additional 
relativistic dark matter species. As a result, the model is viable even for couplings ~ 1. A 
similar mechanism is employed in ref. [11] in order to satisfy the observational constraints. In 
this model the interaction between dark matter and dark energy becomes important only during 
the recent evolution of the Universe. 

It is reasonable to expect that the resolution of the coincidence problem through an inter- 
action between dark matter and dark energy will have to rely on a coupling not significantly 
weaker that gravity. It seems unlikely that a coupling /i^ <C 1 can lead to a cosmological evo- 
lution drastically different from that in the decoupled case. We expect the effects described in 
this letter to be natural consequences of any viable model that achieves the resolution of the 
coincidence problem through the coupling of dark matter to dark energy. 
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